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Abstract 



The optimized linear 5-expansion is applied to the theory at high 
temperature. Using the imaginary time formalism the thermal mass is eval- 
uated perturbatively up to order 5 2 . A variational procedure associated with 
the method generates nonperturbative results which are used to obtain the 
critical temperature for the phase transition. Our results are compared with 
the ones given by propagator dressing methods. 
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I. INTRODUCTION 



The breakdown of perturbation expansion in high temperature quantum field theory is a 
well known problem whose solution is still a matter of study and discussion today, with 
different authors using different methods @-ll|. High temperature perturbation expansion 



breaks down due to the appearance of infrared divergences close to critical temperatures 
(in field theories displaying a second order phase transition or a weakly first order tran- 



sition |HM3|]), or for massless field theories, like QCD. In particular, there are parameter 



regimes where conventional perturbation schemes become unreliable at high temperature 
when powers of the coupling constants become surmounted by powers of the temperature. 

In general, these problems are treated with resummation techniques which try to account, 
in a self-consistent way, for the leading contributions in the infrared region. Among these 
schemes are the popular daisy and super-daisy schemes ||J§], composite operator method 
and field propagator dressing methods [0,0. Some of these resummation methods have been 
compared in Ref. |9|], where their difficulties and possible caveats have also been discussed. 
The majority of approaches used within this subject have a potential drawback concerning 
the achievement of self-consistency as higher order diagrams are resummed. This happens 
because the dressed propagator changes order by order. Therefore, special care must be taken 
when selecting the correct order in the coupling constants. Another problem associated with 
some of these methods is related to the implementation of renormalization, as discussed in 



Ref. flTT 



In this paper, we apply the optimized linear 5-expansion [|l"4| , |T5| (for earlier references see, 
e.g., |L6[]) to the A0 4 theory obtaining the thermal mass to second order in the perturbative 
parameter 5. Our results show that the use of a proper optimization scheme is equivalent 
to solve the gap equation for the thermal mass, where leading and higher order infrared 
regularizing contributions are nonperturbatively taken into account. An advantage of the 
linear ^-expansion is that the same simple propagator is used in the evaluation of any 
diagram, avoiding the potential bookkeeping problems mentioned above. This makes the 
method particularly useful and simple to use in different applications, including the study 
of nonperturbative high temperature effects. 

This work is organized as follows. In Sec. II we briefly describe the linear (^-expansion 
technique and then use it to evaluate the thermal mass up to order 5 2 in the 3+Ig? A</> 4 theory. 
Details of the renormalization up to this order are given in Sec. Ill where we also discuss 
renormalization at higher orders in 5. We show that it does not present any additional 
difficulty when compared to renormalization in the usual perturbative or loop expansions. 
In Sec. IV we present our results for the thermal mass, including the critical temperature 
for the phase transition, and compare them with other results found in the literature. In 
Sec. V concluding remarks are given. 



II. THE LINEAR 5-EXPANSION APPLIED TO THE EVALUATION OF THE 
THERMAL MASS IN THE \<p 4 THEORY 
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A. The Linear ^-Expansion 



The optimized linear 5-expansion is an alternative nonperturbative approximation which 
has been successfully used in a plethora of different problems in particle theory [|T/5|,p~7^[T9 



quantum mechanics pO|j2T|| statistical physics f2~2 |, nuclear matter |23[ and lattice field the- 
ory ||24|| . One advantage of this method is that the selection and evaluation (including 
renormalization) of Feynman diagrams are done exactly as in ordinary perturbation theory 
using a very simple modified propagator which depends on an arbitrary mass parameter. 
Nonperturbative results are then obtained by fixing this parameter. An interesting result 
obtained with this method in the finite temperature domain is given in Ref. |18| where the 
critical temperature value for the Gross-Neveu model in 1+1 dimensions nicely converges, 
order by order, towards the exact result set by Landau's theorem (T c = 0) Q. 

The standard application of the linear 5-expansion to a theory described by a Lagrangian 
density C starts with an interpolation defined by 



C s = (1 - S)C (v) + 5C = C (v) + S[C - £o(v)\, 



(2.1) 



where Cq(j]) is the Lagrangian density of a solvable theory which can contain an arbitrary 
mass parameter (77). The Lagrangian density C 5 interpolates between the solvable Cq{t]) 
(when 5 = 0) and the original C (when 5=1). In this work we consider the A0 4 model 
described by 



2 x 2 j 2 
-mo 



Ct ) 



where 



ct 



A l -{d,cj>f-B l -m 2 cj> 2 



(2.2) 



(2.3) 



represents the counterterms needed to render the model finite. Note that C ct requires an 
extra piece if one attempts to evaluate the thermal effective potential |IT| , which is not the 
case here. Choosing 



and following the general prescription one can write 



C* = ^<P? ~ - 6 ^-i(l- 5)^ + Ci , 



1 

2' 



or 



(2.4) 



(2.5) 



(2.6) 



1 Consistent results in the X(f) theory at finite temperature have also been obtained with a variant 
of the linear (^-expansion, the nonlinear (^-expansion [25,26|. However, beyond first order, this latter 
version presents cumbersome technical problems associated with the evaluation of graphs. 
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where Q 2 = m? + rf and 



C 



ct 



(2.7) 



One should note that the 5-expansion interpolation introduces only "new" quadratic terms 
not altering the renormalizability of the original theory. That is, the counterterms contained 
in C s ct , as well as in the original £ ct , have the same polynomial structure. 

The general way the method works becomes clear by looking at the Feynman rules 
generated by C s . First, the original <p A vertex has its original Feynman rule — iX modified 
to — i5\. This minor modification is just a reminder that one is really expanding in orders 
of the artificial parameter 5. Most importantly, let us look at the modifications implied by 
the addition of the arbitrary quadratic part. The original bare propagator, 



S(k) =i(k 2 -m 2 + ieY 



becomes 



S(k) = t{k 2 - Q 2 + te)- 1 



m 2 + ie 



-i -l 



k 2 — m? + ie 



-17] 



(2.8) 



(2.9) 



indicating that the term proportional to rfcfi 2 contained in Cq is entering the theory in a 
nonperturbative way. On the other hand, the piece proportional to 8rf(j) 2 is only being 
treated perturbatively as a quadratic vertex (of weight i5r] 2 ). Since only an infinite order 
calculation would be able to compensate for the infinite number of {—it] 2 ) insertions con- 



tained in Eq. (|2T9|), one always ends up with a r\ dependence in any quantity calculated to 
finite order in 5. Then, at the end of the calculation one sets 5 = 1 (the value at which 
the original theory is retrieved) and fixes rj with the variational procedure known as the 
Principle of Minimal Sensitivity (PMS) |27| Q 



dP{rf) 
dr] 



0, 



(2.10) 



where P represents a physical quantity calculated perturbatively in powers of S. 

This optimization procedure, together with the convergence problem, has been discussed 
in detail for simple cases in low dimensions in Refs. |21J and []2T| where possible implica- 
tions to more realistic theories have also been investigated. Both references provide proofs 
of convergence. Using the anharmonic oscillator, Bellet et al. 
convergence of an alternative version of the linear (^-expansion. 



2S 



have also studied the 
Their method has been 
later extended to the Gross-Neveu model where the optimization procedure was studied in 
conjunction with the renormalization group [29|. The convergence of the 5-expansion in 



quantum field theory is still a subject deserving further investigation. In principle, it seems 
that if the convergence proof conditions for the anharmonic oscillator could be extended 
to quantum field theory in d < 4 |21], one could pursue a similar investigation for the A0 4 
theory in 3 + 1 d at finite temperature. This is due to the fact that at very high temperatures 



The third refence in [20 discusses alternative conditions for fixing rj. 
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this model gets dimensionally reduced to an effective 3c? theory for the zero Matsubara field 
modes |3(J. However, we shall not pursue this discussion. 

As far as renormalization is concerned it is important to note that in general, as a result 
of the optimization procedure, the arbitrary rj turns out to be a function of the original model 
parameters, scales introduced through regularization as well as external parameters such as 
the temperature and/or density. Therefore, in order to get physically acceptable results the 
optimization procedure must be carried out after all divergences have been eliminated. The 



renormalization problem, in the large- N limit, has been addressed in Ref. |31|. The way 
renormalization will be carried out here is well discussed in Ref. (ITT . 



B. The Evaluation of the Thermal Mass at order 8 2 

We can now start our evaluation of the thermal mass, defined by 

Ml = n 2 + E s T (p) , 



(2.11) 



where is the thermal self-energy. At lowest order (first order in 5) the relevant con- 

tributions, which are momentum independent, are shown in Fig. 1 and given by 



4 1 W = -V + 4/ 7 



d 4 k 



t (2tt) 4 k 2 - n 2 + ie 



(2.12) 



The temperature dependence can be readily obtained by using the imaginary time formalism 
prescription (see, e.g. M) 



p — > lUJ n , 

Then, the self-energy becomes 



It (2tt 



Ef(p) = -V + ^E/ 



rf 3 k 
(2tt) 3 



d 3 k 



1 



(2tt) 3 uj 2 n + E 2 ' 

where E 2 = k 2 + Q 2 . Summing over Matsubara's frequencies one gets 



E?(p) = -6 V 2 + 6- f-^l—- - t 

T w ' 2 J (2vr) 3 \ 2E E[l - ex] 



exp(E/T)] 



Then, using dimensional regularization [32] one obtains the thermal mass 



Ml = tt 2 - 5r] 2 + 5 



A 



32tt 2 



:ti 2 



h In 



ll 2 



+ SXTh ( 5 



where /i is a mass scale introduced by dimensional regularization and 



h(y) 



1 



47T S 



dx 



xr 



[x 2 + y 2 } 2 [exp(x 2 + y 2 ) 2 — 1] 
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(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2-17) 



where x = k/T. Note that the temperature independent term diverges and must be renor- 
malized. In this paper we chose the Minimal Subtraction (MS) scheme where the countert- 
erms eliminate the poles only. At this order the only divergence is 



XQ 2 
32n 2 e 



(2.18) 



which is easily eliminated by the 0(8) mass counterterm 



sS(^ 2 



B sl n 2 



8- 



X 



32vr 2 e 



(2.19) 



By looking at Eq. fl2.16p one can see that the terms proportional to 5X represent exactly the 
same diagram which appears at O(X) in ordinary A0 4 theory, excepted that we now have 
Q 2 instead of m 2 and 5 X instead of A. Therefore, it is not surprising that to this order the 
renormalization procedure implied by the interpolated theory is identical to the procedure 
implied by the original theory at 0(X). 

Let us now analyze the temperature dependent integral which is expressed, in the high 
temperature limit (fi/T <C 1) 



as 



%) = 24-8^-16^ 



V 



"'if 



(2.20) 



In principle, since rj is arbitrary, one could be reluctant in taking the limit Q/T <C 1. 
Therefore, to be sure that the PMS can be safely applied to the thermal mass in the high 
temperature limit, we have performed numerical calculations using both forms for the in- 
tegral h(y) finding that the optimization results do not lead to any significant numerical 
changes. Then, taking the integral h(y) in the high temperature limit, one obtains the 0(8) 
thermal mass: 



M 2 



n 2 -8rf + 8 



XT 2 
~24~ 



8 



XTQ 



8 



in 



XQ 2 
32^2 



In 



'AttT 2 



IE 



(2.21) 



At 0(8 2 ) the self-energy receives contributions from momentum independent as well as 
momentum dependent diagrams. At this order there are five diagrams contributing, which 
are shown in Fig. 2. Let us first consider the momentum independent diagram given by 
the first diagram in Fig. 2, which we call Yi\ . In the high temperature approximation it is 
given by 



Ef 



8 



2 XTrf 



-5< 



32tt 2 




- IE 



(2.22) 



As will be shortly seen, this contribution can be rendered finite using a mass type coun- 
terterm contained in C ct which is tailored to account for divergences arising from the extra 
quadratic vertex introduced during the interpolation process. 

Considering the 0(8) mass counterterm used to eliminate the divergence in E^ iv (see Eq. 
( p.!9| )), one is able to build a one loop 0{5 2 ) diagram whose contribution is given by the 
second diagram in Fig. 2. In the high temperature approximation one obtains 



5 



9 A 2 f2 2 r9 A 2 
-5 2 — _ - + <r 



(327r 2 ) 2 e 2 
A 2 ft 2 



TO fi 2 
+ 



2(32tt 



2^2 



In 



327r 2 e [ 16vr 32vr 2 

n 2 \ 



In 



'47TT 2 



7e 



47T/i 2 



7e 



2 7T 2 

+ 1T 



(2.23) 



Next, one considers the vertex counterterm, whose Feynman rule, — 3i(5 2 A 2 /(327r 2 e), can be 
obtained, as in Ref . [0 , by evaluating the order-5 2 contribution to the four point function. 
The one loop graph evaluated with this counterterm, given by the third diagram in Fig. 2, 
gives 



Ef ~ 



3A 2 f2 2 3A 2 
-5 2 — — - + 5 2 



(32vr 2 ) 2 e 2 



T 2 Ttt tt 2 

+ 



32vr 2 e 24 8tt 32vr 2 



In 



AttT 2 



-JE 



3\ 2 tt 2 



2(32tt 



2^2 



In 



Q 2 



Any 2 



+ 1 



Ti- 



(2.24) 



The next momentum independent contribution is given by the first two loop diagram 
shown in Fig. 2: 



Ef ~ 6 2 



x 2 n 2 



(32tt 2 ) 2 e 2 



A 2 1 T 2 



3TQ 



+ 



n 2 



2 \2 



-<5 2 A 
+ S 2 

+ 6 2 



T 3 



2\2 



384vrfi 
A 2 



32vr 2 e [24 16tt 16tt 2 
T 2 



In 



4ttT 2 



7s 



(16tt) 2 I 

\ 2 n 2 

2Y2 



1287T 2 

T 2 3TQ Q 2 
~3 2t7 + Al 2 



In 



4ttT 2 



(32tt 



In 



n 2 



Any 2 



+ IE 



3 7T 

4 + "6 



IE 
2 



In 



4ttT 2 



-JE 



(2.25) 



To render this diagram finite one needs mass and vertex counterterms 

The final contribution to the self-energy at 0(5 2 ) comes from the two-loop "setting sun" 
diagram shown by the last term in Fig. 2. This is a momentum dependent contribution 
which is given by the real part of 



A 2 

S 2 — (G + Gi + G 2 

D 



(2.26) 



where Go is the zero temperature part (in Euclidean time) of the diagram and G\ and Gi are 
the finite temperature ones (with one and two Bose factors, respectively). Re [Go] is given 
by(d = 4-2e): 



Re[G (p)] = y 



d d k 



d d q 



1 



(2n) d J (2n) d k 2 + Q 2 q 2 + Q? (p - k - q) 2 + 9? ' 
This contribution has been evaluated in details in Ref. MM where the quoted result is: 



2.27) 
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Re[G (p)] 



(27rj" ,: sm rr 



\2d 



2) fe, n = 



2/i 



ttj nlT(\ 



■x 



x 



r(| + fe)s(i + fc, 1 + fe)r(2 - f + fe) r(2 + fe)s(3 -f + M- | + *O r ( 4 - d + fc) 



(Jfe-n)!r(| + ib-n) 



(Jfe-n + l)!r(3 + Jfe-n-f) 



/v 4e 



(27T) 



2c/ 



i ,„r(3 - rf + n) /p\ 2n / d d \ , 



where B(x,y) is the Beta function: B(x,y) = [T(x)T(y)]/T(x + y). In what follows we 
evaluate the self-energy on-shell (p = 0, po = —ifl). For e — > 0, we obtain the following 
result for the above expression: 



A 2 X 2 A 2 \ 2 2 
' ;J 4(4vr) 4 



6 

Wft 2 
"2(4tt) 4 



1 3-2 T£ 2 / Q 2 
-77 H m 



47T/i 2 



4(47r) 4 6e" 



In 



, , Q 2 \ / 17\ , / Q 2 



ATTfi 2 



+ 71-^ + 3.5140 
o 



(2.29) 



where we purposefully left the momentum dependence in the relevant divergent term to 
make explicit the need for a wave-function renormalization counterterm. 

The finite temperature contributions G\ and G2 are given, as in Ref. H, by 



where 



A 2 

- S 2 —Re[G 1 (-in, 0)] = F + F x + F 2 



(4tt) 2 e VT 



(2.30) 



(2.31) 



and 



j2 A 2 T 2 /!) 



In 



ll 2 



47T/i 2 



+ 2- 7s 



d/c 



dq 

E k {eP E »-l)Ja E~ 



gin 



X 4 



1/,- 



(2.32) 



(2.33) 



X± 

For G 2 , one has 
5 2 A 2 



^ 2 — (Ek + E q + Ek± q ) 2 £l 2 — (—Ek + E q + Ek±q) 



(2.34) 



5 2 \ 2 n 2 



6 



Re[G 2 (-ifi, 0)] = H(n) 



dx 



dy 



4(2vr) 4 ii e^ x - 1 ii - 1 

In the high temperature limit, one can show || that 



In 



y/x 2 - 1 + vV - 1 



(2.35) 
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F 2 + H(0) ~ 5' 



A 2 T 2 



24(4vr) 



In 



'ft 2 

y2 



+ 5.0669 



(2.36) 



Finally, using the high temperature approximation for h(y) and putting all together one 



gets 
Repf(p)]^ 2 



A 2 ft 2 1 r2 A 2 ft 2 1 
+ 5 2 



,2 AY 
+ <r 



(32tt 2 ) 2 e 2 



(32tt 



2\2 



A- 



^ Ttt Q 2 
16vr 2 e 1 24 ~ 87 + 32^ 



In 



(32tt 2 ) 2 6e 
/4vrT 2X 



x 2 n 2 



s 2 " 



2(4tt) 4 

2 



+ 5' 



327T 2 

A 2 T 2 
24(4tt)^ 



In 2 
In ( 
In 



47T/i 2 y 

I T 2 



17 



+ (2 7B --)ln 



-7e 

/ n 2 



j + 2 - lE 
+ 5.0669 



v 47T/i 2 

r 2 

24 ~ 87 ~ TeT 2 



ll - ^ + 3.5140 



6 



In 



f— ) 

V4vrT/ 



' IE 



(2.37) 



III. ON THE RENORMALIZATION AT ORDER 5 2 AND AT ORDER ^ 

To obtain the total finite order 5 2 contribution one can add all divergences appearing in 
Eqs. (|2.22|) - (p.25| ) an d (|2.37|) . As it can be easily seem all the nonrenormalizable temperature 
dependent divergences cancel exactly and one is left with 

^S 2 ^S 2 / 2\ , ^<5 2 /„2\ , ^<5 2 



where 



siv = KUW) + KUp 2 ) + KUv 2 ) , (3.i) 



S£v(fi 2 ) = * 2 7^(-^ + ") • ( 3 -2) 



(32vr 2 ) 2 V e 2 e 



S1(P 2 ) = . (3.3) 



and 



<M) = ^\ , (3.4) 

By looking at all diagrams which contribute to this order one can identify two classes. The 
first is composed by diagrams such as the ones described by Eqs. fl2.23j )-( p.26j ). All of them 
are analogous to the diagrams which appear at 0(X 2 ) in the original theory and can be 
rendered finite by similar mass and wave-function counterterms, which are respectively 

SS(^) = B*'tf = (-1 + i) I? , (3.5) 
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and 



The second kind of diagram is exclusive of the interpolated theory and carries at least one 
5r] 2 vertex. At 0(5 2 ) this diagram is described by Eq. ( p.22| ) which displays the divergent 
term Y*f iv (r} 2 ). Looking at C 5 ct one identifies a rj 2 counterterm whose Feynman rule is i5B 5 rj 2 . 
Since the actual pole is of order 5 2 one then identifies the coefficient as being B 5 , displayed 
in Eq. ( gTg) . Then, 

nt(v 2 ) = -^ 1 V 2 = -6 2 -^ V 2 . (3.7) 

In practice, the renormalization at higher orders can be done as above. That is, 0(5 n ) 
diagrams belonging to the first class will be renormalized exactly as in the original theory 
at 0(\ n ). This is obvious from the fact that all the diagrams in this class are of order- 
S n X n . It is easy to check that for those diagrams the most divergent terms will display e _n 
poles. On the other hand, 0(5 n ) diagrams belonging to the second class will make use of 
the counterterm SB S " r/ 2 , where B s n has been evaluated in a previous order. One can 
also easily check that for these diagrams the most divergent terms will have e n_1 poles. 
Moreover, power counting reveals that those St] 2 insertions make the loops more convergent. 
For example, all one loop diagrams of order 0(5 n ), with n > 3 are finite. 

Finally, one should note that the renormalization prescription adopted here is in accor- 
dance with the one suggested in Ref. [IT|, where the order by order renormalization was 
shown to hold at any higher orders in 5. 



IV. NUMERICAL RESULTS 

One can now set 5 = 1 and apply the PMS to the finite thermal mass. First let us set 
m = so that our results for the thermal mass can be compared directly with the resummed 
perturbative expansion (RPE) results of Ref. [[| . At order-5 one gets 

(4-1) 

which, clearly, does not depend on the the coupling and cannot generate nonperturbative 
information. However, nonperturbative results appear already at second order in 5. Table 
1 shows our results and the results furnished by the RPE for M|. in units of /i for A = 0.1. 
We also show, in units of the optimal values of 77. 



77 = 2nT 
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T/im 


T}/fi 


0(5') 


RPE 


0.5 


0.033 


0.098 


0.099 


1.0 


0.067 


0.391 


0.396 


1.5 


0.100 


0.880 


0.892 


2.0 


0.133 


1.564 


1.587 


2.5 


0.166 


2.445 


2.481 


3.0 


0.200 


3.522 


3.574 


3.5 


0.233 


4.794 


4.867 


4.0 


0.266 


6.263 


6.358 


4.5 


0.299 


7.927 


8.049 


5.0 


0.332 


9.788 


9.939 


5.5 


0.366 


11.845 


12.029 


6.0 


0.399 


14.098 


14.317 


6.5 


0.432 


16.547 


16.806 


7.0 


0.465 


19.193 


19.494 



Table 1 - Results for M 2 //i 2 (xlO 



-2\ 



Let us now obtain the critical temperature for the phase transition at 0(S 2 ). Taking 
A = 0.1, we reset m 2 = — /i 2 in M|> observing a second order phase transition at the critical 
temperature T c = 15.57 /i whereas the modified perturbation scheme (MPS) of Banerjee and 
Mallik [0] predicts T c = 15.63 /i. Choosing A = 0.01 we find T c = 49.03 \x whereas the value 
T c = 49.05 fi is predicted by the MPS. Note that in the calculation of the critical temperature 
performed in Ref. |7j the propagator has been effectively dressed up to the leading order 
correction in the temperature, which is set by the tadpole term in fl2.12| ) (see their Eq. 
(4.5)). Here, on the other hand, we are definitely working with all higher order corrections 
up to the two-loop contribution. The fact that our results for the critical temperature are 
slightly smaller than those obtained in Ref. is an indication of the importance of these 
higher order corrections and is in accordance with well known results concerning the study 
of phase transitions in the context of the electroweak effective potential beyond 1-loop ||13| . 
The results are also in accordance with recent results for the finite temperature effective 
potential of the A0 4 theory, obtained with the super-daisy approximation [55]. 



V. CONCLUSIONS 

Using the A0 4 model we have shown how the optimized 5-expansion can be useful in 
extracting nonperturbative information through an essentially perturbative evaluation of 
Feynman graphs. Our 0(5 2 ) results for the thermal mass and for the phase transition critical 
temperature are in excellent agreement with the ones given by other methods |7]|| . However, 
although providing very similar results, these methods differ, from the ^-expansion, in some 
aspects which may become important if one tries to consider higher orders. For example, 
within the latter methods the effective mass used in the modified propagator changes order 
by order turning the propagator into a coupling dependent quantity from the start. One 
can then expect the selection and evaluation of higher contributions to become complicated 
quickly. The ^-expansion method avoids these potential problems by using Q in the modified 
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propagator which is used at any order calculation. Therefore, after drawing the relevant 
graphs which contribute to a given order, one does not have to worry about bookkeeping 
inconsistencies nor renormalization problems since this is done as in perturbation theory. 
The extension to higher order in S is immediate and, as discussed above, leads to a consistent 
resummation procedure in finite temperature field theory. 

Although we have not attempted to prove the possibility of convergence of our results we 
have explicitly shown that the procedure interpolation, renormalization, and optimization 
in the finite temperature domain can be consistently handled to furnish encouraging results. 

We also note that the linear 5-expansion can be extended to the case of gauge theories, 
where it has already been used as a tool to study the electroweak phase transition on the 
lattice |p4|| . Recently, it has been shown |36] that the method does not spoil gauge invariance. 
In this context, the linear ^-expansion may be a useful technique to analytically study the 
nonperturbative aspects and difficulties associated, for example, with the electroweak phase 
transition as well as other problems in high temperature gauge theories. 
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Figure Captions 

Figure 1: Diagrams contributing to the self-energy at first order in 5. 
Figure 2: Diagrams which are order 5 2 contributing to the self-energy. 
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FIGURES 
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